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Abstract 

We show how to obtain the first multipole contributions to the electromagnetic radiation emited 
by an arbitrary localized source directly from the Jefimenko equation for the magnetic field and the 
Panofsky-Phillips equation for the electric field. This procedure avoids the unnecessary calculation 
of the electromagnetic potentials. 
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I. INTRODUCTION 



The derivation found in most textbooks of the electromagnetic fields generated by arbi- 
trary sources in vacuum starts by calculating the corresponding electromagnetic potentials 



(see for instance, Refs. , ory). After the retarded potentials are obtained (assuming 
the Lorenz gauge) the electromagnetic fields are calculated with the aid of the relations 
E = — V*!* — {1/ c){dA/ dt) and B = V x A (we use Gaussian units). The resulting ex- 
pressions for the fields are usually called Jefimenko's equations because they appeared for 
the first time in the textbook by Jefimenko.- Jefimenko's equations are obtained in Ref. Il 
from the retarded potentials and are obtained directly from Maxwell's equations in Ref. |5. 
(Heras^ had already derived Jefimenko's equations directly from Maxwell's equations.) Ref- 
erences y andH obtain a less common form of Jefimenko's equations for the electric field, 
but this form is more convenient for studying radiation. 

Griffiths and Heald- illustrate Jefimenko's equations by obtaining the standard Lienard- 
Wiechert fields for a point charge. Ton^ provides an alternative derivation of Jefimenko's 
equations and three applications, including that of a point charge in arbitrary motion. 
Heras has generalized Jefimenko's equations to include magnetic monopoles and obtained 
the electric and magnetic fields of a particle with both electric and magnetic charge in 
arbitrary motion. He has also discussed Jefimenko's equations in material media to obtain 
the electric and magnetic fields of a dipole in arbritary motion^^ and has derived Jefimenko's 
equations from Maxwell's equations using the retarded Green function of the wave equation.- 

The main purpose of this paper is to enlarge the list of problems that are solved directly 
from Jefimenko's equations (or the equivalent). Our procedure avoids completely the use 
of electromagnetic potentials. Specifically, we shall obtain the electric dipole, the magnetic 
dipole, and the electric quadrupole terms of the multipole expansion due to the radiation 
fields of an arbitrary localized source. 

II. JEFIMENKO'S EQUATIONS FROM MAXWELL'S EQUATIONS 

In this section we present three methods of calculating Jefimenko's equations directly from 
Maxwell's equations. The first method closely follows Ref. 0. The second method makes 



use of a Fourier transformation as discussed by Ref. [l6|. For our purposes it suffices to do 
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a Fourier transformation only in the temporal coordinate. We shall see that this method, 
although longer than the previous one, avoids any possibility of misleading manipulations 
with retarded quantities. The subtleties in calculations involving retarded quantities have 
been discussed in Refs. Iisljljjisl . We then present an alternative method that is a variation 
of the first one; the main difference is the order of integration. 



A. Direct calculation using the retarded Green function 

The following approach is similar to that in Refs. 0, and 13. Maxwell's equations with 
sources in vacuum are given by 

V ■ E = 47rp (1) 

V ■ B = (2) 

VxE^-if (3) 

c ot 



Taking the curl of Eq. ([3]) and the time derivative of Eq. (Ill), we obtain 
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)E.4.(V,+ If). (5) 

where we used Vx (VxF) = V(V-F)- V^F, with V^F = {V'^F.^)^+{V'^Fy)y+{V'^F,)i 
and F being an arbitrary function. Similarly, we obtain for the magnetic field: 

V=-iS)B^-i^X (6, 



The solutions of Eqs. ([SD and ^ can be obtained with the aid of the retarded Green function 
Greti^,t;x.' ,t'), which satisfies the inhomogeneous differential equation 

1 

and is zero for t — t' < 0. The solution for Gret for t — t' > is given by^ 

r (^f^'f^- 1 ^(t--(t-|x-x-|/c)) 1 6{t'-{t-R/c)) 

where R = |R| = |x — x'|. With the help of Eq. ([8]) the solution of Eq. ([5]) can be written as 

1 aj(x',t') 



E(x,t) = -J^J dt'5{t' - (t - R/c)) (^V'p(x',t') + 



c2 dt' 



(9) 
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where the notation [. . .] means that the quantity inside the brackets is a function of x' and 



is evaluated at the retarded time t' = t 



ic. 



At this point much care must be taken, because V'[p] ^ [Vp]. Due to a bad choice of 



notation, the author in Ref. 



13| used, incorrectly, the quantity [V'p] as if it were V'[p] and. 



after an integration by parts, an incomplete result for the electric field was found, as pointed 



out in Ref. 

Ref, 
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The correct relation is given by V'[p] = [V'p] + R-[p]/c (see, for instance. 



. [s]). If we use the correct relation, Eq. ffTOl) becomes 



E(x,t) 



rfx'M^ + 



(hi' 



cR 



dx.' 



c^R 



dx 



cR 



dx.' 



c^R 



(11) 
(12) 



In the last step we integrated by parts and discarded surface terms, because the charge 
distribution is localized. Equation (fT2!) is one of the Jefimenko's equations.-i^ Note the 
retarded character of the electric field. The first term on the right-hand side of Eq. f|T2l) is 
the retarded Coulomb term. 

As shown by Panofsky and Phillips,- there is an equivalent way of deriving Eq. fll2p for 
the electric field which manifestly shows the transverse character of the radiation field. A 
discussion of this point can be found in Ref. [s]; we will discuss this point in Sec. III. 

To obtain the desired expression for the magnetic field, we use the retarded Green function 
in Eq. ([8]) to rewrite Eq. ([6]) as 



Bfx 



(13) 



R 



If we use the relation V' x [J] = [V' x J] + — x [J] (see Ref. Q), Eq. ([T3]) takes the form 



B,x,0 4/=|V'X|J|-J,/. 

We integrate by parts and obtain 



,R X [j] 
R~ 



B(x,*) 4/<;x' V X + i/.x' Ul X V'(i) - i/ 



dx' 



,R X [j] 



R 



(14) 



(15) 



The first term on the right-hand side of Eq. f|T5|) vanishes because the current distribution 
is localized in space. We use the relation V'(l/-R) = — R/i?^ to obtain 



B(x,t) 



dx 



[J] X R [J] X R 

ci?2 + c^R 



(16) 
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The first term on the right-hand side of Eq. (fT5|) is the retarded Biot-Savart term; the 
transverse radiation field is contained in the last term. Equations f|T2l) and (|T6|) are the 
Jefimenko equations. 



B. Fourier method 

In Sec. Ill Al we showed how to obtain Jefimenko's equations by a careful treatment of 
the derivatives of retarded quantities. This point is crucial - spatial derivatives cannot 
be commuted with retarding the functions, because the retarded function depends on the 
coordinates in its time argument. A simple way to circumvent this difficulty is to use Fourier 
transforms and factor out the time dependence of the functions so these subtleties are not 
encountered. We will show how to obtain the electric field and will leave the derivation of 



the magnetic field as an exercise for the interested reader (see, for example, Ref. Il6l). 

We start with the electric field given by Eq. (Q. The term involving the current density 
takes the form after integration over time. 

We introduce the Fourier transformation p(x',co') of p(x',t') as p(x',t') = / p(x', ct;)e~*'^* du, 
and express the remaining contribution to the electric field in Eq. (ITOi) as 

-j^J dt'Kt' -it- R/c))v'p{^',t') 

= -J^J dt'6{t' - (t - R/c)) J rfcuV'p(x',cu)e-*"*' (18a) 

dooe-'^' / dx'^ V'p(x', uj) (18b) 

= -jdue-^-'j dx'|v'(^p(x',.;))-p(x',..)(^-zfc|)e^^^| (18c) 

= j rfx'-^ j rfcup(x', cu)e^'"(*-^/^) + j d^'^J cia;p(x',o;)(-^cu)e-'"(*"^/^) (18d) 

To obtain Eq. f ll8cp we integrated by parts; to obtain Eq. ( ]18d[) we discarded the surface 
term. 

We combine Eq. with Eqs. (fT7|) and (fT8|) and obtain the electric field generated by 
arbitrary (but localized) sources, namely, Eq. (fT2|) . On the right-hand side of Eq. (I18ap it is 



obvious that V' acts only on p{x.',uj). Hence, in Eq. f llScP there is no possibihty of thinking 
that V' also acts on the exponential e*^^., The subtleties of dealing with retarded quantities 
are circumvented by the Fourier method. 



C. Postponing the delta- function integration 

When we are faced with integrals involving delta functions, we are usually tempted to do 
them first, because they are so easy. To derive Jefimenko's equation for the electric field, 
given by Eq. f[T^ . this procedure is not optimum. Let us start with Eq. (jH]) 

E(x, t) = -j^j dt'5{t' - (t - R/c)) ( V'p(x', t') + lg(x', t')) . (19) 

Instead of first performing the time integration using the Dirac delta function, we integrate 
by parts on both terms on the right-hand side of Eq. f|T9|) so that it takes the form 

E(x,t) = j dx'j dt'l^6'{t' - (t - R/c)) (^i-p(x',t')V'i? + -^J(x',t')) 

+ 5(t'-(t-i?/c))(v'i)p(x',t')|. 
If we use the relations V'i? = — R and V'l/i? = R/i?^, we obtain 

E(x,t) = I dx'l dt'i^6'{t' -it- R/c)) |^--^p(x',t') + ;^J(x',t': 

+ 5(t'-(t-i?/c))p(x',t');||- (21) 

By using the properties of the Dirac delta function, we can easily perform the integration 
over t' to obtain Eq. f[T^ . An analogous calculation provides an expression for the magnetic 
field. 



(20) 



III. MULTIPOLE RADIATION VIA JEFIMENKO'S EQUATIONS 

The main purpose of this section is to add to the list of problems that can be handled 
directly with Jefimenko's equations by calculating the first multipole contributions to the 
radiation fields of an arbitrary localized source. We shall obtain the first three contributions. 
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namely, the electric dipole, the magnetic dipole and the electric quadrupole terms. Most 
textbooks treat this problem by calculating first the electromagnetic potentials. 

Although Eq. (IT^ for the electric field gives the correct expression for the electric field 
of moving charges, it is preferable for our purposes to write it in an equivalent form as given 



in Refs. 



7 and 



E(x,t) 



+ / rfx'M:^)^ + ([J] X R) X R 



.x'MiiS^. (22) 



i?2 7 cR^ J c^R 

We obtain Eq. ( l22l) starting with Jefimenko's equation for the electric field in Eq. ( fT2ll . 

Sl. Note that (the Einstein convention of implicit 



Our derivation will follow the one in Ref. 



summation over repeated indices is assumed) 



V ■ f Jl 



[ V • J] + [I] 



d 



t'=t-R/c 

cR 



■Ip] + (j) 



t'=t-R/c 

R 

c 



(23a) 
(23b) 



where Xi = Xi — x[ and in the last step we used the continuity equation. From Eq. fl23bp 
we have [p] = —V' ■ [J] + [j] • R/c, which can be substituted into the second term on the 
right-hand side of Eq. ([T2|) to yield 

'v'-iJllR 



[p]R 



dx' 



-dx' 



1 



J R R 



-dx'. 



(24) 



cj R cj R c^J R 

We show that the first term on the right-hand side of Eq. ( 12^ is proportional to 1/R'^, so 
that it does not contribute to the radiation field. Observe that 



R 



C 

■^fd'J[j.]§)d.'+^l[m 



Xk 



dx.' 



2 fJl - R R 



R^ 

Jl R^R 



-dx.' 



i Jl X R X R 



i?2 



-dx\ 



(25a) 
(25b) 
(25c) 

(25d) 

(25e) 



where we have used the identity (a x b) x c = (a ■ c) b — (b ■ c) a, and the fact that the 
surface term appearing in Eq. (125bp vanishes because the sources are localized. We substitute 



Eq. (125|) into Eq. and insert the result into Eq. (fT2l) to obtain 
E(x. *) = / 4^ dx' + i/ (|Jl-R)R+mxR)xR 

.x-M + /,x-'W ■ ^) + /dx-'W f ^. (26b) 

R'^ / cit^ / c^R 



which is Eq. (1221) . If we consider arbitrary time- varying sources at rest and use Eqs. (11611 
and (|22|) . the (transverse) magnetic and electric radiation fields are given by 

B,ad(x,t) = - / d^' ^ , Erad(x,t) = / d^^^^^ • (27) 

It can be shown^"^ that for time- varying sources in motion, Eq. (!27|) gives the radiation fields 
plus additional non-radiative terms of order 0{1/R'^). For instance, if a Hertz dipole is 
accelerated, integration of equations Eq. ( |271) yields radiation fields plus nonradiative terms 
of order C(l/i?^) which are induced by the dipole motion (for details, see Refs. llTI ). 

In the radiation zone we can write R ~ x, 1/R ~ 1/r and i? ~ r — x • x', where we 
defined r = |x|. If we substitute these approximations into Eq. (1271) . we obtain 

B,ad(x, t)c^— d^'j (^x', to + X X (28) 

X /* X x' 

Erad(x, ^) - ^ / f^x' [ j (^X', to + X x] X X, (29) 



where = t — r / c is the retarded time of the origin. A comparison of Eqs. (!28l) and (j3j 
leads to the relation 

E,ad(x,t) =Brad(x,t) X X. (30) 

Now we are ready to calculate the first multipole contributions for the radiation fields. 
We need to calculate only one of the radiation fields because the other is readily obtained by 
Eq. fl30l) . which also shows that the fields in the radiation zone are mutually orthogonal. We 
start by calculating the electric dipole term and then consider the next order contribution 
given by both the magnetic dipole and electric quadrupole terms. 



A. The electric dipole contribution 



The lowest order contribution to the radiation fields comes from the electric dipole term. 
For simplicity, we calculate the lowest order contribution to the radiation magnetic field, 
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which we denote by B^^^, 



B«(x,t) = ^^ /rfx'j(x',to)[^xx. (31) 



c^r 



We write the unit vectors of the cartesian basis as e, = V'x^, (i = 1,2,3), and write any 
vector V as V = e^Vi = ej(v ■ e^). The integral in Eq. flHT]) can be expressed as 

j rfx' j(x', to) = Gij d-x! j(x', to) ■ = Gij rfx' j(x', to) ■ V'x^ (32a) 

= rfx'V ■ {x'i j(x', to)) - e^j dx.' x[ V ■ j(x', to), (32b) 

where in the last step we integrated by parts. Because we are considering localized sources, 
the first integral on the right-hand side of Eq. f l32ap vanishes (after the use of Gauss' theorem 
this integral is converted to a zero surface term). The remaining integral may be cast into 
a convenient form if we use the relation 

V'.J(x'.t„) = -^!4^. (33) 

which is a direct consequence of the continuity equation. To obtain Eq. ( I33l) we used the 
relation 9p(x',to)/9t = 9p(x', to)/(9to (see Ref. Il8|, note 5). We substitute Eqs. ( |33l) and 
fl32al) into Eq. (13T]) and obtain 

p(*o) 

B:ikx,i) = i^{/ rfx'p(x'.i„)x'}xx=EM^, (34) 

where p(to) is the electric dipole moment of the distribution at the retarded time to. Now 
it is clear why this first term is called the electric dipole term. The radiation electric field is 
readily obtained from Eq. (!30|l 



Eax.t)= '— (35) 

Equations ( l3Ti) and ( l35i) are the radiation fields of the electric dipole term, that is, the 
first-order contribution to the multipole expansion. These expressions are valid for arbitrary 
but localized, sources in vacuum such as an oscillating electric dipole. 

B. Next order contribution 

To calculate the next order term we need to take into account the second term of the 
expansion 

j(x',to + ^^) j(x',to) + ^^J(x',to). (36) 



[p(to) X x] X X 



We substitute Eq. (136|1 into Eq. (128!) and identify the next order contribution to the 



(2) 

radiation magnetic field, o^J^, given by 

BS(x,t) „3 



1 

c^r 



rfx' J(x',to)(x ■ x) I X X. 



(37) 



For reasons that will become clear we will split the integral in Eq. fl37|l into antisymmetric 
and symmetric contributions under the exchange of J and x'. This rearrangement will give 
rise to the magnetic dipole and electric quadrupole terms of the multipole expansion for the 
radiation fields. 

If we write J(x', Iq) (x ■ x') as - J(x', to) (x • x') + 2^(-^'' ^o)(x • x') and sum and subtract 
-[J(x', to) ■ x] x' in the integrand of the right-hand side of Eq. ( 1371) . we obtain 



B;;i(x,t) 



c^r 2c 



+ 



dx' 



J(x',to)(x ■ x) - ( J(x',to) ■ x)x' 



X X 



c^r 2c 



dx.' 



J(x',to)(x-x') 



J(x',tn) -x x' 



X X. 



(38) 



For pedagogical reasons we shall treat the magnetic dipole and electric quadrupole cases 
separately. 

Consider the antisymmetric term of Eq. ( l38l) . We denote this contribution by B^^^'*, use 
(a X b) X c = (c ■ a)b — (c ■ b)a, and write it in the following suggestive way 



B 



{2a) 
rad ' 



1 1 



c^r 2c 



c^r 



x,t) = ^ / fix' J(x',i(:o)x ■ x' - J(x', to) ■ X X' 



X X 



fix' 



,x' X J(x',to) 



2c 



X X > X X, 



(39a) 
(39b) 



m(io) 



where we have identified the second time derivative of the magnetic dipole moment of the 
charge distribution at the retarded time as m(to). Hence, the antisymmetric contribution 
for the radiation magnetic field is given by 

[m(to) X x] X X 



B2?(x,t) 



(40) 



The appearance of the magnetic dipole moment of the distribution justifies the name given 
to this contribution. The corresponding radiation electric field is readily given by 

X X m(to) 



Efa?(x,t) 



c^r 



(41) 
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Expressions (HOj) and (HT!) are valid for an arbitrary, but localized, time-varying source at 
rest in vacuum such as an oscillating magnetic dipole. 

We denote the symmetric term of Eq. (155]) as B^^^ . We have 



b2?(x, t) = j;Ycl '^'^ ^ ^^"^^ ■ ""'^ + ( ^^"^^ to) • X )x' ^ X X. 



(42) 



We manipulate the first integral on the right-hand side of Eq. fj42l) in the same way as before. 
We have 



j rfx'J(x', to) (x ■ x') = Bij rfx' (x • x') J(x', to) ■ 



ei /(ix' (x-x')J(x',to) ■ Vx^. 



If we integrate by parts and remember that the surface term vanishes, we obtain 



j rfx'J(x',to)(x ■ x') = -e,y" rfx'x^V ■ (x ■ x')J(x',to) 



(43a) 
(43b) 

(44a) 



V'(x-x')) • J(x',to) + (x-x')V'- J(x',to) . (44b) 



-Cj / (ix' x' 



We use V ■ J(x',to) = -d^p{x.',to)/dt^, the relationi^ dp{x.' ,to)/dt = (9p(x', to)/<9to), and 
V'(x ■ x') = X to obtain 

fix' J(x',to)(x ■ x') = -Gj /"(ix'x ■ J(x',to)x- + /"c/x' (x ■ x')a;-p(x', to), (45) 



which implies 



J(x', to)(x ■ x') + ( J(x', to) ■ X )x' = ^ / rfx'(x ■ x')x'p(x', to). (46) 



The left-hand side of Eq. f l45l) is the integral in Eq. fH2]) . We write the factor 1/2 in Eq 
as 3/6 and obtain 

1 r 



B2?(x,t) = i rfx'(3x. x')xV(x',to)| X X 



6c^r 
1 I 



6c^r I dt^ 



(47a) 
(47b) 



where in the last step we included in the integrand the term r'^ xp(x', to), where r' = |r'|, 
which gives a vanishing contribution to the result because x x x = 0. Now, we define the 
transformation Q such that 



Q(£,t)=y"rfx'[3(|-x')x' + r'^|] p(x',t). 



(48) 
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This transformation takes a vector and an instant of time {$,t), and transforms it into the 
vector Q(^, t), given by Eq. fHHj) . If we use the definition (j^Hj) . Eq. fl47al) for the antymmetric 

(2) 

contribution for becomes 

B2?(x,t) = ^Q(x,to)xx. (49) 

The corresponding symmetric contribution for the radiating electric field is obtained from 
Eq. ([30]): 

Eiad (x, ^) = ^ [Q ^o) X x] X X. (50) 

In order to interpret the above term, let us define a linear operator Q*, for a fixed instant t, 
by Q*(^) = Q(^, t). Note that Q* is a linear operator so that Q*(ai^i + 02^2) = <^iQ*(^i) + 
a2Q*(^2)5 '^15 G The linear operator Q* is called the electric quadrupole operator. 
Since Q*(ej) is a vector in 7^'^, we can write it as a linear combination of the basis vectors, 
namely 

3 

Q*(e,) = 5^Q;,e, (^ = 1,2,3). (51) 

i=i 

The coefficients Q*j are the cartesian elements of the electric quadrupole tensor (a second 
rank tensor) of the distribution at instant t. That is why we interpret Eqs. (H9l) and (l50!l as 
the electric quadrupole contribution for the radiation fields. 

The electric quadrupole radiation fields given by Eqs. ( l49l) and (!50l) . together with the 
magnetic dipole radiation fields given by Eqs. f HOl) and ( HTl) . are the first corrections to the 
leading order term given by Eqs. and (1351) . In this sense, we can write the first multipole 
contributions to the multipole expansion for the radiation fields of a completely arbitrary, 
but localized, time-varying source at rest in vacuum as 

Brad(x,t) = ^ |p(to) X x + [m(to) X x] X X + ^ Q (x, to) X X + . . .| (52) 

Erad(x, t) = Brad(x, t) X X. (53) 

With these radiation fields, we can calculate the corresponding Poynting vector and the 
power radiated by the source. In this paper we have presented an alternative approach 
of obtaining the first multipole contributions to the radiation fields of general radiating 
systems directly from Maxwell's equations. We hope that it may be useful for those who 
enjoy studying electromagnetism. 
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